Vibrational modes and spectrum of oscillators on a scale-free network 
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We study vibrational modes and spectrum of a model system of atoms and springs on a scale- 
free network in order to understand the nature of excitations with many degrees of freedom on the 
scale-free network. We assume that the atoms and springs are distributed as nodes and links of a 
scale-free network, assigning the mass Mi and the specific oscillation frequency uii of the i-th atom 
and the spring constant Kij between the i-th and j-th atoms. 
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I. INTRODUCTION 



Recently there has been a notable progress in the study 
of the so-called scale- free network (SFN) 1, 2, 3,,.4, 5, 6, 7, 
Hi 1^ ^1 ■ In the point of view of network theory, the ran- 
dom network theory (RN) was first invented by Erdos 
and Renyi^ and has been applied to many areas of 
sciences fro m p hysics such as Anderson localization[T^ . 
percolation ^J, and energy landscape 0, ^| to biology 
such as the KaufFman's NK modelflal- Recently it was 

Eeralizcd to the small world network (SWN) models 
J3^M2d^^^2£2M and has been applied 
to physical systems such as percolation 0, 0, 0[ I20L 
I21II22II23II23 125^ . arrays of coupled oscillators with syn- 
chronization known as the Kuramoto model.26>] , arra ys o f 
coupled lasers l27l, arrays of the Josephson junctions [23, 
electric circuit '29l, traffic transportationj^^ |^, and 
protein folding 32J . And very recently, the SFN was 
discovered by st udy ing the network geometry of the 
internetlE B IS IssL l34 Issl Is^ . Albert, Jeong and 
BarabasilzH. 133 . 134 ISSllSq opened up an area for study- 
ing very coinplex and growing network systems such as 
internet^, |33|, .34 , 35, .36,. .,37.J . biological evolution ,1 3^. 
metabolic react ionl39l . epidemic disease [ioj. human sex- 
ual relationship|4l|. and economy These are nicely 
summarized in the reviews by BarabasiQ- 

The nature of these SFNs is characterized by the 
power-law behavior of the distribution function for the 
number of nodes with k links such as P{k) oc fc~^ where 
7 « 1 — 4. In order to show the power-law distribu- 
tion of the SFN, Albert and Barabasi first proposed a 
very simple model called the Albert-Barabasi (AB)'s SFN 
modeljl H m HI m 113. This system is constructed 
by the following process: Initially we put mo nodes as 
seeds for the system. Every time when a new node is 
added, m new links are distributed from the node to the 
existed nodes in the system with a preferential attach- 
ment probability Ili{ki) = ki/ ^f^^^ ki, where ki is the 
number of links at the i-th node and we have assumed 
m < toq. The development of this model is described 



by a continuum model ^ — mlli{ki) — Then at 

time r the system consists of N{t) nodes and the L{t) 
links with L(r) = ^ 'J2illi'' ^i- This model exhibits 7 = 3 



for the power-law. Thus, it has been concluded that the 
essential points of why a network grows to a SFN are at- 
tributed to the growth of the system and the preferential 
attachment of new nodes to old nodes existed already in 
the network. 

However, although time evolution of the SFN has 
been intensively studied regarding nodes and links as 
metaphysical objects such as agents and relationships 
in an area of science, it seems that very few physi- 
cal models putting real meaning on nodes and links 
in the SFNs have been studied in order to investi- 
gate excitations such as vibrations, phonons, and elec- 
trons, except diffusion and spins 13 on the SFN 
and excitations in the RN 11, 12, iXlll HIL and 
SWN[i1 El E I23, 0, m m . So, we explore 
to study, as a prototype model, vibrational modes and 
spectrum of a system of atoms coupled by springs where 
the atoms and springs are located regarding as nodes and 
links of a SFN. 



II. VIBRATIONAL MODEL 

Let us introduce our vibrational model. We first adopt 
the AB's SFN model for the construction of a SFN and we 
regard nodes and links in the SFN as atoms and springs 
in our physical model. Assuming that qi and Wi are the 
displacement and the specific frequency of the i-th atom 
of mass Mi, respectively, we can define the hamiltonian 
of the system: 



H 



E 



M, 



M,ujf 



Nir) 

E 



(1) 

where — ^ the velocity of the i-th atom and Kij is 
the spring constant between the i-th and j-th atoms with 
Kij = Kji. Although we will concentrate to study only 
this atom-spring model in this paper, the generalizations 
of this model and the applications to other systems are 
straightforward. We expect that the physical nature of 
this model shares with those of such models and systems. 

We now assume that the time interval Ar for the de- 
velopment of link addition process is much larger than 



2 



that At of the physical model such that |At| ^ |At|. 
This guarantees that although the network grows in the 
course of its development, as long as the network consists 
of N{t) nodes and L{t) links, the vibrational model can 
be simultaneously solved. This means that time evolu- 
tion the network is adiabatic to the time motion of the 
atoms and springs. By using the Euler-Lagrange equa- 
tioii. ^(ff) = ff> we obtain 

N{t) 

for i = 1, . . . , N{t). Assuming q^{t) = q,{uj)e-'-'^\ Eq.(2) 
becomes 

N(t) 

M,iu;f - u;^)q, = ^ A>(q, - q,), (3) 

for i = 1, . . . ,N{t). This is the eigenequation for our 
system. 

Let us assume that all springs are identical for the 
sake of simplicity such that Kij = KgAij, where Kq is 
the spring constant and Aij is the ij-th component of the 
adjacency matrix A for the network geometry. The com- 
ponents of the adjacency matrix are non- negative such 
that Aij = or 1 according to whether or not a link be- 
tween the i-th and j-th nodes exist in the network. The 
link number ki{T) at the i-th atom (i.e., the order of the 
j-th node) is given by kiir) = X^jL^i"* ^ji- From this, the 

last term in Eq.(3) becomes -^jili — Koki{T)qi. 

Hence, in this setting, we obtain 

N{t) 

fl,q, = KoJ2 ^J^lj^ (4) 
i=i 

for i = 1, . . . , N{t), where 

n, = M,{ujf-u;^)+KohiT). (5) 



= Go + i^oGoAG = Go + GoTGo, (7) 
where T is called the T-matrix defined as 

t = i^oA + A'o^AGoA + • • • = ti + ta + • • • . (8) 
We can now derive the following: 

(t„),, =ifo"(AGoA•••GoA),, 
~ ^ o- o. ...o. — = J^ohj . (yj 

and since Kij — Kji (i.e., Aij = Aji), we find 

(t„),, - (t„),,. (10) 

From Eqs.(8) and (9), we find (t)y = KoVl]'' + K^T^f + 

Let us consider the trace of the Green's function. We 
now get 

rrG = Tr(Go + GotGo)=^(^^ + ^^. (11) 

We note here that if Mi{ujf — w^) = such that = 
Koki{T), then (T„)ij = KohPij with fc^Py = kjPj,, 
where Pij (r) means the probability that the walker starts 
at node i at time t = and found at node j at time t = t 
in terms of the language of the diffusion theory of Noh 
and Rieger [see Eqs(2) and (3) in j4^]. As is well-known, 
the density of states p{uj) is given by 

p(cj) = --TrG(cj + ie). (12) 

Thus, the poles of the Green's function produce the spec- 
trum of the system. 

IV. SPECIAL LIMITS 



III. GREEN'S FUNCTION FORMALISM 

Let us now define the Green's function by 

N{t) 

[^iSij - KoAij] Gjk = Sik, (6) 

for i, k — 1, . . . , N{t), which is represented by [Gq ^ — 
KqA]G = 1 in the matrix representation where 1 is the 
N{t) X N{t) unit matrix and Go is the N{t) x N{t) diag- 
onal matrix defined by Gq = ^~^Sij. Thus, the Green's 
fimction is formally obtained as G^^ = Gq"^ — ii'oA. Fur- 
thermore, we can derive a series expansion of G in terms 
of Go and A as 

G = Go + /^oGoAGo + i^o'GoAGoAGo + • • • 



Before going to do the direct calculation for spectrum 
of the system, let us consider some limits, (i) The in- 
dependent atom limit. First, in the case of no springs 
of Kq = 0, since Tr(G) = Tr (Go), the poles of the 
Green's function are given by ^lj = Mj{uXj — uP') — 0, 
which trivially provides the discrete spectrum uj ~ ujj for 
j — 1,...,7V(t). This means that the atoms indepen- 
dently vibrate with specific frequencies uji. 

(ii) The AB limit. Second, in the case of very weak 
spring constant such as Kq ^ 1, the the poles of the 
Green's function are obtained as 

+ (13) 

This means that each atom vibrates with frequency re- 
lated to the number of links of the atom. Since the distri- 
bution of the nodes with k links is given by P{k) oc k~^ 



in the SFN0,0], the distribution of the spectrum is given 



as 



(14) 



Hence, this hmit shares with the nature of the AB's SFN 
geometry. Therefore, we may call this limit the AB limit. 

(iii) The localized mode limit. Third, let us consider the 
limit of very small mass (Mi ^ 1) or very strong spring 
{Kq ^1). In this case, we can ignore the frequency 
dependence in the eigenequation of Eq. (4) such as fli ~ 
Koki, which then yields 



9j> 



(15) 



for i — 1, . . . , N{t). Since we can rewrite the above equa- 
tion as '^j'J:i\kiSij — Aij]qj — 0, non-trivial solutions 
may exist only when the determinant det[ki6ij — Aij] 
vanishes. This is realized when qj = qi where j runs 
the adjacent links around the i-th atom. In this sense 
the mode is localized within the adjacent atoms. 



V. CALCULATIONS OF THE SPECTRUM 

Let us now calculate the spectrum of the system of 
oscillators in the SFN. This is carried out by directly di- 
agonalizing Eq. (4) . For the sake of simplicity, we assume 
that Mi = Mq — const, and uji = ujq = const, and we 
adopt the AB-model for generating the SFN geometry. 
We have performed the calculations for the systems up 
to = 10^. 

Fig.l (a) shows the density of states of the system, 
where we have calculated for the case of to = 2 and 
N = 10^ (blue) and the case of m = 4 and Af = 10^ 
(red), respectively, where we have used toq = 5. To ob- 
tain the distributions, we have used twenty configurations 
with different random numbers. The vertical axis means 
pi^)\/{k)2, while the horizontal axis means oj'^/ ^ {k)2, 
where (fc)2 stands for the second order average degree of 
a node [See Eq.(18)]. The shape of the curve is unique 
such that there is a peak at w = cjo and the spectral tail 
exists in the whole range of the spectrum. This tendency 
means that there is a scale-free nature in the spectrum 
of the vibrational modes in the system. 

Fig.l (b) shows the tail behavior of the density of 
states. The density of states is shown in a log-log plot 
for the cases of the AB-model with to = 2 and A^ = 10^ 
(blue circles), with to = 4 and N = 7 x lO'^ (red crosses), 
and with to = 4 and A^ = 10^ (red squares), respectively. 
The red line is a guide for showing (uj'^)~^. 

From this we find that the tail behavior of Eq.(14) 
holds valid for the general cases as well. Therefore, we 
can conclude that our vibrational model shares common 
nature with the AB-model of the SFN. This is contrary 
to the conclusion previously obtained from the calcula- 
tions of the spectrum of the adjacency matrix A of the 
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FIG. 1: (color online) The density of states and its tail behav- 
ior of the system of oscillators in the SFN. (a) The density of 
states is calculated for the cases of the AB-model with m — 2 
(blue) and m = 4 (red) for N = 10*, respectively. The distri- 
butions are obtained as an average over twenty configurations 
with different random numbers. The vertical axis means the 
density of states, p{uP)^ {k)2, while the horizontal axis means 
oj^, where u is the vibrational frequency of the oscillators and 
{k)2 the second order average degree of a node, (b) The tail 
behavior of the density of states is shown for the cases of the 
AB-model with m — 2 and N — 10* (blue circles) , with m = 4 
and N = 7 X 10^ (red crosses), and with m = 4 and N = 10* 
(red squares) , respectively. The vertical axis means log-plot of 
the density of states while the horizontal axis means log-plot 
of uj^, where lu is the vibrational frequency of the oscillators. 
The red line is a guide for showing {uj'^)~^. Here we have 
assumed that loq = Ko = Mo = 1 and mo = 5. 



AB-model |45l l46l l47l l4a | . There, when the network has 
the tail behavior of P{k) (x k~'^, the spectral tail for 
the eigenvalues A of the adjacency matrix is given by 
p(A) oc A^'' where 7' = 27— 1. Therefore, since the AB- 
model has 7 = 3, we conclude 7' = 5. This is different 
from our result of 7' = 3. The main reason for this phe- 
nomenon is explained as follows: In our vibrational model 
the fii consists of the degree ki of the node [see Eq. (4)] . 
Therefore, as the system grows, so does the magnitude of 
fli. This can reduce the contributions of the adjacency 
matrix A in the higher terms of the perturbation series 



4 



of Eq.(7). Hence, the spectral behavior is dominated by 
the pole of the unperturbed Green's function Gq. Thus, 
we are led to the same spectral behavior in the AB-limit. 

The physical meaning of the above results can be un- 
derstood as follows: The main peak in the density of 
states is attributed to vibrational modes with frequency 
ujQ. These modes are extremely localized within the least 
connected nodes in the SFN such that the total number 
of the localized modes provides the height of the peak. 
Since the number of modes is nothing but the number of 
degeneracy of the eigenequation, these localized modes 
are highly degenerate. 

On the other hand, there is the power-law tail of 
p{u)'^) oc {uj"^)^^ as u! oo. This means that the larger 
the frequency of modes the fewer the number of modes. 
In other words, as the frequency is increasing, the num- 
ber of modes is decreasing by the power-law. As the re- 
sult, there appears only one mode with the maximum fre- 
quency (i.e., the maximum eigenvalue). The mode with 
the maximum frequency is extended over the entire sys- 
tem of the SFN. This situation means that in the SFN 
the lowest frequency modes can be very easily excited, 
but it is very hard to excite the maximum energy mode. 
Thus, the high frequency modes are very hard to exist 
in the system of oscillators coupled in the SFN. This na- 
ture is very different from that of the standard systems 
of networks such as rn|ii1 and lattices 13 that there 
are a small number of orders of nodes. This is the most 
prominent characteristic of our system. 



VI. THE MAXIMUM EIGENVALUE 

The behavior of the maximum eigenvalue Xmax of the 
adjacency matrix A is very important in the network 
theory 1, 2, 45l l46Ll47l l4l L In the standard networks such 
as the random networks ITU , the maximum eigenvalue 
Amoij: cannot grow so fast as the network grows [3, 0, EE 
|43,|43,|43. And also, as in solid state physics, networks 
in most of physical systems provide the so-called energy 
band that is a spectrum with a finite region [T3.iT^ . This 
is due to the topology of the finite coordination number of 
atoms in the network of the lattice structure [l^ . So, 
in order to elucidate the difference between the SFNs and 
other networks the growth of the maximum eigenvalue is 
an important signature. 

As was numerically studied by many authorslH B S 
I45I IS . 23 1 the maximum eigenvalue Xmax of the ad- 
jacency matrix A in the AB-model is proportional to 
^/kmax such that 



A. 



cx 



(16) 



Here kmax means the maximum order of nodes in the 
network such that kmax ~ maxijfci} (We will use this 
notation for later purposes). And the numerical studies 
showed that kmax \/N . Therefore, we obtain 



To see whether or not this is true in an arbitrary SFN 
and to know how general it is, very recently, Chung, Lu 
and Vu48!] have proved a very general theorem: 

Theorem 1 Suppose that the distribution of degrees of 
nodes in a SFN is represented by P{k) cx k~^ . Denote 
by (fc)2 the second order average degree of a node. This 
is defined by 



N(r) 



E 



N{t) 



(k) 



(18) 



where (k^) — ^y^Ty EiL''i''' kf with p integer. Then, (CI) 
if the exponent 7 > 2.5, then 



const. {k) 2 < Xmax < const.y/hn 



(C2) If the exponent 2 < 7 < 2.5, then 



const. Vk 

max < Xmax < const, (k) 2. 



(19) 



(20) 



( C3) And if the exponent 7 = 2.5, then a transition hap- 
pens. 

We note here that in the paper of Chung, Lu and Vu|43| 
they used the notation d for the second order average 
degree, instead of our notation (fc)2 for it. Applying the 
above theorem to the AB-model of 7 = 3, we find that 
the AB-model belongs to the first category. Hence, the 
theorem explains the numerical results [a. Isl lisl l46l l47| . 

In spite of such efforts, whether or not the growth of 
the maximum eigenvalue of a physical model on the SFN 
is not so well-known. This is because the eigenvalues 
of the adjacency matrix is different from those of the 
eigenequation of a physical system. In this sense, the 
problem to investigate the growth behavior of the maxi- 
mum eigenvalue of the eigenequation of a physical system 
is a nontrivial problem. So, in order to see this point, let 
us consider the maximum eigenvalue LOmax of our vibra- 
tional system of oscillators. 

We have performed calculations of the maximum eigen- 
value (i.e., vibrational mode) Wmax in our model of oscil- 
lators on the AB-SFN, where m = A and N is developed 
up to = 10"*. This is shown in Fig. 2. The maximum 
eigenvalue cj^^^ (circles) , the maximum degree of a node 
kmax (triangles), and the second order average degree 
(fc)2 of a node (-I-) are shown, respectively. Here we have 
obtained the following relation: 

iol + 2{k)2 < < u;l + 2kmax. (21) 

This looks similar to the result of Eq.(19) such that 



(22) 



A^ 



cx 



(17) 



However, this is not supported by our numerical calcu- 
lations. Therefore, as the spectral tail of our vibrational 
model is different from that of the AB-model as discussed 
in the previous section, so is the growth behavior of the 
maximum eigenvalue of our vibrational model. This is 
an important character of our physical model with the 
AB-SFN. 
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FIG. 2: (color online) The behavior of the maximum eigen- 
value in the spectrum. The maximum eigenvalue cj^ax (cir- 
cles), the maximum degree of nodes luq + 2{k}2 (triangles), 
and the second order average degree ujq + 2kmax of nodes 
(squares) are shown in the vertical axis, respectively. The 
horizontal axis is scaled as m^/iV(r). The calculations have 
been carried out for (a) the AB-model with m = 4 and A'^ 
is up to Ai" = 10*, and for (b) the AB-model with m = 2 
and A^ is up to A = 10*. Here we have taken the values of 
Ko = Mo = uoQ = 1 and (fc)2 stands for the second order 
average degree of a node. 



VII. SOME THEOREMS ON THE MAXIMUM 
EIGENVALUE 



Let us consider the origin of the inequality, Eq.(21). To 
see this, let us go back to Eq.(4). From the Hadamard- 
Gerschgorin's theoremfl^ (see also Appendix A) we can 
derive an inequality 



N{r) 



(23) 



Since \qj\/ 



\qi\ < 1 and \Aji \ = Aji, we can derive |r2i| < 
I = Koki{T), which then yields a theorem: 



Theorem 2 



2 2 



Thus, there exists at least one atomic site (i.e., node) 
that satisfies Eq.(24) for all eigenvalues uj. This implies 
that llP' is included within a disk of radius ^fci(T) and 

its center iJ^ -I- jf-ki{T). 

Since in our model all Kij (i.e., Aij) are non-negative, 

(see also 

-cofl < 
satis- 



by applying the Perron-Frobenius's theorem|lr 
Appendix B) to Eq.(24) we can derive |a;^ 
2^/ci(T). Hence, the maximum frequency w^^^ 
lies another theorem: 



Theorem 3 



< max 



(25) 



And similarly we can obtain a more precise theorem: 
Theorem 4 



< ^max < max 



(26) 



From this, if we assume that Mi = Mq — const, and 
uji = ujQ — const, and applying for the SEN, then Eq.(25) 
becomes 



Theorem 5 

,^0 



2^min{fc,(r)} < c^^.-^o < 2-5^ max{A;,(T)}. (27) 

Hence, this theorem verifies our numerical results in the 
previous section. Therefore, the upper limit of the spec- 
trum (i.e., spectral edge) grows as fast as the network 
grows. This is a remarkable fact for excitations in the 
SEN models and this nature is very different from that 
of Anderson localization where only mobility edge may 
appear in the spectrum and the band edge cannot grow 
as fast as the system size grows [T^ IT^ . 

The above Theorems 2-5 are good for the standard 
networks that the distribution of the orders of nodes is 
limited such as periodic lattice systems or the RN[13>| or 
the SWM[ll[ll,[li|23,|2ll[22, 23, 24111, since in these 
systems there exist finite lower and upper limits of the 
orders of nodes such that the error width is bounded as 



A(a;: 



2 

max 



-^o) > max[2^A:j(r) 



-min[2^fe(r)]. (28) 

I Mi 



However, whether or not the above theorems can be ac- 
curate conditions for the SEN 0,13 is not trivial, since in 
the SENs there exist various orders of nodes without any 
bound but with the power-law distribution. 

To study this point, we first observe that there is a 
particularly important nature of the adjacency matrix A 



6 



in the network theory. Denote by ki = {An, . . . , ^iAr(7-))* 
the i-th column vector of A. The vector represents the 
way of hnks between the i-th node and other Unked 
nodes, such that it defines the order ki of the i-th node 
such that 



ki = kl ■ ki — y ^ Ai 



(29) 



Therefore, let us call ki vectors the link vectors. Using 
this representation, we can rewrite the adjacency matrix 
as A = {kl,. . .,kN(r)) = {ki, 
the transpose. From this, we can derive that 



fc^j-^P*, where t means 



(fcj • kj). 



(30) 



which is a symmetric matrix and nothing but the 
Gramian matrix between the link vectors, fc, , where 



Therefore, for the maximum frequency we obtain 



uif I < max 



A/. 



k,(T) 



Hence, we obtain 



^'Ln-r < ^i + max 



Since k. 
max. 



< 



k* 



max. 



^(fc,(r) + yfc^ 
f (^fc,(r) + ^4,.fc,) 



(37) 



(37) 



ki , we have 

< 

Therefore, the right 



hand of Eq.(37) is comparable with that of Eq.(26). In 
this way, Theorems 2-5 work for the SFN systems as 
well. 



N{r) 



N{t) 



TT{K')=^kl-h^ ^ 

1=1 1=1 



h - 2L{t) 



(31) 



Let us go back to Eq.(4). We now rewrite it as ^iq — 
KoAq. Therefore, fifg = K^A?q. Let us now use the 
Hadamard-Gerschgorin theorem[l3 (Appendix A) or the 
Perron-Frobenius theorem^^ (Appendix B) for KqA^, 
we can derive an inequality 



\nj\<K, 



Nir) 
2 \ ^ I 



Since < 1 and [(A^ 

derive m<K^j:A{A'\ 



(32) 



we can 



ki, where 



'^tot 



Then we have 



Nir) 

= ^fc* = (fci,fc: 



■ , kN)- 



2l.t 
tot 



i=l 



< 



M, 



(33) 



(34) 



Therefore, it then yields a theorem: 
Theorem 6 



Ko 



h{T) 



< 



Ko 



ktot ■ ki 



(35) 



Thus, there exists at least one atomic site (i.e., node) that 
satisfies Eq.(35) for all eigenvalues w. This implies that 

is included within a disk of radius ^ \/ ^tot ' ki and 



its center ujf + ^ki{T). Since \uj'^ 



f'-ki{T] 



< 



^h{r) 



we obtain 



< ^ ( kAr) + 



(36) 



VIII. CONCLUSIONS 

In conclusion, we have studied the system of oscilla- 
tors connected by springs in the geometry of the AB SFN 
model. We first presented the Green function formalism 
for obtaining the spectrum of the vibrational modes of 
the system. In the case of very weak spring constant, us- 
ing this formalism we find that the distribution of eigen- 
modes follows the same type of power-law distribution 
of degrees of a node in the AB model [see Eq.(13)] such 
that P{uj'^ — Wq) oc (tj^ — ojo)~^ with 7 = 3. In the case 
of an arbitrary strength of spring constants, we have per- 
formed numerical calculations in order to obtain the spec- 
trum of vibrational modes. We have found that even in 
this case, the distribution of eigenmodes obeys the same 
type of the power-law distribution of degrees of a node 
in the AB model as well [see Fig.l]. This is contrary 
to the distribution of eigenvalues of adjacency matrix in 
the AB model, where power-law distribution is given by 
p{\) (X A"''' with 7' = 27 - 1 = 5. 

This is a consequence of our model, where relative 
displacements between the individual oscillators are in- 
cluded in the Hamiltonian. This Hamiltonian provides 
the diagonal matrix elements in the eigenequation, which 
are proportional to the degrees of nodes [see Eq.(4)]. 
These diagonal elements can be regarded as on-site po- 
tentials in the problem. Since the degree of a node de- 
velops indefinitely, the on-site potential can be arbitrary 
large as the system is progressing. Therefore, the eigen- 
values are strongly dominated by the magnitude of the 
diagonal elements of the eigenequation. Thus, the dis- 
tribution of eigenmodes is affected by that of degrees of 
nodes such that the distribution of eigenmodes coincides 
with that of degrees of nodes in the network. 

We finally have investigated the asymptotic behavior 
of the maximum eigenvalue ujmax of the system. We 
have found numerically that the maximum eigenvalue is 
bounded as in Eq.(21). From this, as the total num- 
ber of nodes, N, is increasing, the maximum degree of 
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nodes becomes arbitrarily large. Therefore, the maxi- 
mum eigenvalue can be arbitrarily large as TV — > oo. This 
coincides with the result of the maximum eigenvalue of 
adjacency matrix in the AB model. We have also proved 
the above numerical results by some mathematical the- 
orems that are proved using the Hadamard-Gerschgorin 
theorem and the Perron-Frobenius theorem. 

Thus we conclude that when we apply a certain phys- 
ical model to the geometry of a SFN, the physical prop- 
erties are strongly dominated by the nature of the SFN. 
In this sense, not only the network geometry of a SFN 
but also the property of physical models on a SFN are 
important in the study of the SFN. This direction will be 
very interesting for further researches. 
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for i = 1, . . . , iV. Then, we find 

\\-h,\< f2 \h,M, (A2) 

for i = 1, . . . , A^. Since always < 1, we obtain 

AT 

\X-h,\< ^ |/iy|=B„ (A3) 

for i — 1, . . . , N. Now, we find a theorem that there exists 
at least one site such that the above equation Eq.(A.3) is 
valid for all A. Eq.(A.3) means that A is included within 
a disk of radius Bi with its center of hi . 



This paper is dedicated to the memory of Dr. Mihoko 
Yoshida (Lehigh University) who always helped us for 
collecting relevant papers. We would like to thank Dr. 
Jun Hidaka for sending us relevant papers. K. I. would 
like to thank Kazuko Iguchi for her financial support and 
encouragement . 

APPENDIX A: THE 
HADAMARD-GERSCHGORIN THEOREM 

The following theorem is known as the Hadamard- 
Gerschgorin theorem in linear algebra p^. Consider the 
following eigenequation: 

TV 

{X-hi)q,= ^ hijQj, (A.l) 
3^H¥^ 



APPENDIX B: THE PERRON-FROBENIUS 
THEOREM 

The following theorem is known as the Perron- 
Frobenius theorem in linear algebra 0. Suppose that 
an n X n symmetric matrix H has all non-negative en- 
tries hij > 0. Then this satisfies an eigenequation 
H\ijji) = Xilipi). For any positive constants ci, C2, . . . , c„, 
the maximum eigenvalue Xmax{H) satisfies 

>ymax{H) < max \ V I . 
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